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Interactions between cold ions and atoms have been proposed for use in implementing quan- 
tum gates 1 , probing quantum gases 2 , observing novel charge-transport dynamics^, and sym- 
pathetically cooling atomic and molecular systems which cannot be laser coolecP^. Fur- 
thermore, the chemistry between cold ions and atoms is foundational to issues in modern 
astrophysics, including the formation of stars, planets, and interstellar clouds^, the diffuse 
interstellar bands 7 , and the post-recombination epoch of the early universe 8 . However, as 
pointed out in refs 9 and 10, both experimental data and a theoretical description of the ion- 
atom interaction at low temperatures, reached in these modern atomic physics experiments 
and the interstellar environment, are still largely missing. Here we observe a chemical reac- 
tion between ultracold 174 Yb + ions and 40 Ca atoms held in a hybrid trap. We measure, and 
theoretically reproduce, a chemical reaction rate constant of K = (2 ± 1.3) x 10 10 cm 3 s 1 
for 1 mK < T < 10 K, four orders of magnitude higher than reported for other heteronu- 
clear cases. We also offer a possible explanation for the apparent contradiction between 
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typical theoretical predictions and measurements of the radiative association process in this 
and other systems. 

Of particular interest in cold ion-atom physics are binary charge exchange (CEX) reactions, 
which limit the performance of the aforementioned proposals and are important inputs into the 
astrophysical models. Binary CEX reactions occur via three distinct mechanisms: (i) non-radiative 
charge transfer (nRCT), A + + B ->■ A + B + , (ii) radiative charge transfer (RCT), A + + B ->■ 
A + B + + 7, and (iii) radiative association (RA), A + + B — > (AB) + + 7. Due to the scarcity of 
reaction data and full quantum calculations, a number of approximation techniques for these three 
processes have been developed, e.g. approximating all RA rate constants by 10~ 14 cm 3 s _1 (ref. 11), 
Landau-Zener theory 12 , the Demkov coupling method-^, and the semi-classical optical potential 
method-^. These techniques, which have varying degrees of accuracy depending on the dominant 
CEX pathway, must be experimentally verified as demonstrated by a recent measurement where a 
factor of ~ 100 discrepancy was found with experiment ! 10 * 15 ! 

Total CEX rate constants have been reported for two ultracold ion-atom species^S. In the 
first system the CEX rate constant for the resonant, homonuclear case of Yb + + Yb was mea- 
sured by Grier et al. to be ~6xl0~ 10 cm 3 s~ 1 (ref. 16). On the other hand, a measurement of 
the rate constant for the non-resonant, heteronuclear case of Yb + + Rb by Zipkes et al. yielded 
~3.5 x 10~ 14 cm 3 s -1 (ref. 17); an experiment on Ba + + Rb appears to have yielded similar results^. 
These values agree well with conventional wisdom for resonant versus non-resonant CEX. Inter- 
estingly, although the RA process is typically predicted to be the dominant mechanism for low- 
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temperature, non-resonant CEXpH^l, none of these experiments observed the formation of a molec- 
ular ion. Further, Zipkes et al. found the ratio of the rate of nRCT to RCT to be K nRC T^BCT ~ 
2.3, while most calculations predict K nRC T^RCT < 10" 3 in the ultracold regimePEH. 

Here we report a measurement of the CEX chemical reaction rate constant, K, for the non- 
resonant, heteronuclear case of 174 Yb + + 40 Ca in a hybrid trap which is four orders of magnitude 
higher than reported for other heteronuclear cases. In addition, by monitoring the production of 
CaYb + molecular ions in the hybrid trap, we place an apparent upper bound on the branching ratio 
of the RA process. Molecular potential curves for the CaYb + molecule have also been constructed 
and used in a calculation that explains this enhanced rate constant as the result of an avoided 
crossing and strong transition dipole moment. In what follows, we describe the experimental 
setup, results, and theoretical model. We conclude with a possible explanation for the apparent 
contradiction between typical theoretical predictions and measurements of the RA process in this 
and other experiments^. 

Figure 1 shows the hybrid magneto-optical and ion trap (MOTION trap) system used in 
these experiments 5 . Yb + ions, created by laser ablation of a solid Yb target, are trapped in a linear 
quadrupole trap (LQT)^and laser-cooled using laser beams (A = 369 nm, 935 nm) aligned along 
the trap axis. In the same region, ultracold Ca atoms are produced and held in a magneto-optical 
trap (MOT) constructed with laser beams (A = 422 nm, 672 nm) that intersect at the LQT center 
(See Methods). Fluorescence images of both the Ca MOT and the ion cloud from two separate, 
nearly orthogonal camera angles allow a 3D reconstruction of the ion and atom clouds and the 
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Figure 1: MOTION trap schematic. Yb + ions and Ca atoms are laser cooled and trapped in the 
central region. Typical images of the ultracold clouds, shown in the upper left of the figure, are 
used for a 3D reconstruction that yields the degree of overlap. Ions are also detected with a voltage- 
gated CEM. 
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degree of overlap is quantified as $ = / pc a (r)pi(f)df , where p Ca is the unit-peak normalized 
Ca atom density and pj is the unit-integral normalized Yb + density (See methods). An example of 
such images and the reconstruction is shown in the inset of Fig. 1 . 

In absence of the Ca MOT, the trap lifetime of the laser-cooled Yb + ions, measured by 
monitoring the ion fluorescence, is ~ 240 s. However, when an overlapped MOT is introduced, 
the ion lifetime decreases to ~ 10 s, as shown in Fig. 2a. This induced ion loss is also confirmed 
visually by imaging the ion cloud during the decay, as well as observing a decreased number of 
detected ions using the CEM. The inset in Fig. 2a shows that the loss of ions is not due to ion 
heating, as such an effect would increase the ion cloud size. From data such as these we calculate 
a loss rate constant using K = (av) = where T is the measured fluorescence decay rate and p 
is the peak density of the MOT. 

Given the large measured rate constant and the disparity between previous work and typical 
theoretical predictions, we performed several experiments to identify possible systematic errors. 
First, in all previous ultracold experiments the rate constant was measured, as it was here, by 
observing the decay of the parent ion, rather than the rate of appearance of the reaction products, 
which are often untrappable. Thus, we set out to verify that no kinematic effects were responsible 
for the loss of Yb + ions from the trap. To this end, we measured the loss rate constant while varying 
only q, the Mathieu parameter, which determines the ions' stability and micromotion amplitude in 
the LQT^, as well as when only varying the trap depth, T D . These data are shown in Fig. 2b. As 
the rate constant appears independent of these fundamental trapping parameters, it is unlikely that 
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Figure 2: a Yb + ion trap lifetime in the presence of the MOT (solid curve) compared to the lifetime 
without a MOT (dashed curve). Dotted curves are single exponential decay fits. Inset: Images of 
the Yb + ion cloud at various times, b Loss rate constant as a function of the Mathieu stability 
parameter q, the trap depth T D , and the excited state fraction of Ca, p ee . 
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the loss is due to a kinematic effect in the LQT. 

In addition, we performed an experiment where BaCl + ions, produced by laser ablation of an 
additional target (see Fig. 1), were sympathetically cooled by co-trapped, laser cooled Yb + ions. 
We then compared, after 20 s of interaction with the MOT, the number of ions remaining in the 
trap and found that while the Yb + ions had decayed away, the BaCl + ions remained. As BaCl + 
ions are energetically forbidden to charge exchange with Ca atomP, we conclude the loss of Yb + 
ions is due to a CEX reaction. 

Finally, as discussed in ref. 23, it is expected that the observed reaction is between ground 
state atoms and ions, since the long-range ion-neutral interaction and ultracold temperatures pre- 
clude excited state atoms and ions from being in close enough proximity to undergo a chemical 
reaction. To verify this, the loss rate constant was measured as a function of the excited state frac- 
tion, p ee , of Ca atoms in the MOT, by varying the cooling laser intensities. As seen in Fig. 2b, 
no dependence was observed and we conclude that the loss of Yb + fluorescence is due to a CEX 
chemical reaction of ground state Ca atoms and Yb + ions. 

After these systematic checks the loss rate constant was measured as a function of temper- 
ature, as shown in Fig. 3. For these measurements the MOT temperature was held at 4 mK and 
the ion temperature varied by changing the cooling laser detuning and intensity. The horizontal 
error bars account for the possible contribution to the average ion energy that results from residual 
radial micromotion^. The observed rate constant is roughly four orders of magnitude larger than 
typical non-resonant, heteronuclear CEX values and is constant within experimental error over all 
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Figure 3: Total experimental (points) and theoretical (black line) charge exchange rate constants. 
Also shown are the calculated rate constants for the three pathways. 



8 



measured temperatures. 



The CEX reaction products can in principle be trapped, but due to the large mass ratio of 
Ca + (m = 40 amu) and Yb + (m = 174 amu) ions, Ca + ions are heated by the cold Yb + ions^. This 
effect is exacerbated by the fact that trap operation in a mutually stable region necessarily puts one 
or both ion species near a boundary of their LQT stability region. As a result, at the trap parameters 
used in this study, Ca + ions have a measured trap lifetime of only ~1 s. Given the comparatively 
low production rate from CEX reactions with Yb + ions, the equilibrium number of Ca + ions in the 
trap is < 2, which is below our current detection limit. Conversely, the mass of the CaYb + reaction 
product is much closer to that of Yb + and does not suffer from these effects. Therefore, neglecting 
further chemical reaction with the Ca MOT, CaYb + ions are expected to have a trap lifetime of 
at least that of uncooled BaCl + ions (m ~ 174 amu) which we measure to be ~30 s. Thus, by 
allowing the Yb + ion fluorescence to decay to background levels and then measuring the number 
of ions remaining in the trap using the CEM, we place an apparent upper limit on the branching 
ratio for RA at 0.02, under the assumption that all detected ions are CaYb + . 

Also shown in Fig. 3 are the results of our theoretical model for the CEX chemical reaction 
rate constant. Since the structure of CaYb + was not previously known, we first calculated the 
electronic potential surfaces that dissociate to the Ca+Yb + and Ca + +Yb limits, as well as their 
transition dipole moment using the MOLCAS software suited, as shown in Fig. 4a and b (solid 
curves), respectively. This approach worked well for the structure of the BaCl + molecular iorPI 

Using these potentials, we calculated the rate constant for the CEX pathways via a quantum 
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Figure 4: a The X 2 S + and A 2 S + potential energy curves and b the corresponding transition 
dipole moment of the CaYb + molecular ion as a function of the internuclear separation R. The 
dotted curves are the diabatized potentials constructed from these potentials. The dashed curve 
gives the coupling constant between the diabatic potentials. The vertical dotted line shows the 
inner turning point at zero collision energy of the A 2 S + potential and the corresponding transition 
dipole moment. 
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mechanical scattering calculation that assumed the particles collide on the A'S + potential and 
charge exchange to bound or continuum states of the X 2 S + potential. Since the nRCT pathway 
occurs via a nonadiabatic transition near the avoided crossing, we determine its rate constant by a 
coupled-channels method with diabatized potentials - black dotted lines in Fig. 4a (See Methods). 
The thermalized nRCT rate constants are shown as a function of temperature in Fig. 3. We have 
varied the strength of the coupling between the diabatic potentials by 20% to characterize the 
nature of the coupling. Although we find that the nRCT rate coefficient dramatically increases 
with increasing coupling strength, it does not exceed a few times 10~ 14 cm 3 s _1 , indicating that the 
scattering is almost adiabatic. 

For the calculation of the RCT and RA rate constants, we use the perturbative Fermi's golden 
rule for the spontaneous emission of a photon by the diatom (See Methods). Emission to the 
continuum of the X 2 S + potential is characterized as RCT, while emission to a bound state of the 
potential is characterized as RA. Figure 3 shows the thermalized radiative rate constants due to 
both RCT and RA mechanisms versus temperature; in both cases on the order of a hundred partial 
waves are included to ensure numerical convergence. We find rate coefficients that are many times 
1CT 11 cm 3 /s and also observe that CEX should produce nearly equal amounts of Ca + atomic ions 
and CaYb + molecular ions. 

While the total calculated rate constant, also shown in Fig. 3, exhibits satisfactory agreement 
with the total measured rate constant, the predicted RA branching ratio of ~ 0.5 strongly disagrees 
with the measured value < 0.02. These ratios could be reconciled, however, if the CaYb + further 
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reacted with the Ca MOT atoms via subsequent RA reactions, e.g. CaYb + + Ca — > Ca 2 Yb + + 7, 
to form heavier molecular ions, which would be unstable in our LQT. Given that RA will pre- 
dominantly produce CaYb + in high-lying vibrational states, which are typically very reactive, this 
explanation seems quite plausible. Further, if the same explanation is applied to the data of ref. 
1, the interpretation of trap loss as evidence of nRCT could be reinterpreted as evidence of RA, 
yielding a RA branching ratio that agrees well with typical theoretical predictions. 

In conclusion, we have presented an experimental and theoretical investigation of the CEX 
reaction of the Ca + Yb + system. The total measured and calculated rate constants are in good 
agreement and are surprisingly large for a non-resonant, heteronuclear CEX reaction. We have 
also offered a plausible explanation for the apparent contradiction between typical theoretical pre- 
dictions and measured rates of RA in this and other experiments. Given the importance of these 
types of reactions for determining astrophysical processes, this work highlights the need for fully 
quantum calculations of atom-ion systems and a renewed effort in laboratory astrophysics. Finally, 
we note that although the CEX rate constant of this system is large, it should not interfere with most 
of the proposed applications of hybrid atom-ion systems. 

Methods 

MOTION Trap. The MOTION trap (Fig. 1) consists of a radio-frequency linear quadrupole trap 
(LQT) ^ (250 kHz < f < 400 kHz, V RF « 300 V) with four sets of segmenting electrodes to 
provide axial confinement. A channel electron multiplier (CEM) mounted above the trap provides 
a species non-specific method of ion detection, and is voltage-gated so that it does not affect the 
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ion dynamics before detection. Yb + ions, created by laser ablation of a solid Yb target with a 
focused ~1 mJ, 10 ns, 1064 nm laser pulse, are trapped in the middle region of the trap and laser- 
cooled using two laser beams aligned along the trap axis; one with a wavelength and intensity of 
A = 369 nm and I fa 400 mW/cm 2 , respectively, to excite the cooling transitions, and one with 
A = 935 nm and I fa 2000 mW/cm 2 to regain ions which fall out of the cooling cycle. 

The ultracold Ca atoms are produced and held in a magneto-optical trap (MOT) using 6 
primary trapping laser beams (A = 422 nm, I fa 4 mW/cm 2 ), a deceleration beam (A = 422 nm, 
7^31 mW/cm 2 ), and a repumping beam (A = 672 nm, / fa 70 mW/cm 2 ), all of which intersect 
at the LQT center. The source of Ca is a chromate-free getter, ~5 cm from the MOT center. The 
typical MOT atom number, density, and temperature are measured by absorption and fluorescence 
imaging and found to be N Ca = 3.7±0.5 x 10 6 atoms, p Ca = 7± 1 x 10 9 cm" 3 , and T Ca = 4± 1 mK, 
respectively. 

Ion number and temperature are determined by observing ion fluorescence while scanning 
the 369 nm laser (/ = 500 MHz/s) and fitting the resultant truncated Voigt profile^. As the ion 
temperature is known, this allows for a determination of the total number of ions in the trap giving, 
among other things, a method to calibrate the CEM. The laser is then slowly scanned up to a 
specific, predetermined red-detuned frequency and is locked via a transfer cavity to the 422 nm 
laser that is referenced to the Ca 4 1 P«— 4 : S transition. 

When operating the LQT in a region where both Ca + and Yb + ions are stable, a large heat 
load on the laser cooled Yb + ions is observed due to the ionization of Ca atoms in the 4 X P state 
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by the Yb + cooling laser. To eliminate this effect, we modulate the intensity of the 422 nm and 
369 nm lasers out of phase. Additionally, this effect is eliminated if the LQT is operated where 
Ca + ions are not trapped. 

Alignment of MOT and Ion Clouds and Overlap Calculation. The position of the ions can 
be shifted axially in the trap by applying an asymmetric bias to the DC segmenting electrodes, 
and the location of the Ca MOT center can be changed both by adjusting the magnetic field and 
relative laser intensities. These controls allow the ion and atom cloud positions to be precisely 
overlapped, which is the starting point for the data reported here. We quantify this overlap, in 
order to determine the effective density (p$) of neutral Ca, by imaging both clouds with two, nearly 
orthogonal cameras. The images, along with the known distances and angles of the cameras, allow 
a full 3D projection of cloud density distributions in the lab reference frame. We then determine 
the overlap factor by convolving the unit-peak normalized MOT density distribution with the unit- 
integral normalized ion cloud distribution. In addition, we time average the overlap factor by 
taking into account the diminishing ion cloud size in coincidence with the fluorescence decay (see 
Fig. 2a). 

nRCT Calculation. We determine the nRCT rate constant by a coupled-channels method with 
diabatized potentials. For the coupling between these potentials, we assume a Lorentzian in R 
with a maximum value equal to half the minimum splitting of the adiabatic potentials. The width 
of the coupling matrix element is chosen such that diagonalization of the 2x2 potential matrix 
reproduces the original adiabatic potentials. Due to the strength of the atom-ion interaction, up to 
several hundred partial waves are required for numerical convergence of the nRCT rate constants 
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at our temperatures, which leads to a large number of shape resonances in the unthermalized rate 
constants as a function of collision energy. 

RCT and RA Calculation. For the calculation of the RCT and RA rate constants, we use the 
perturbative Fermi's golden rule for the spontaneous emission of a photon by the diatom. In prin- 
ciple, both the initial and final state wavefunctions are solutions of a two-channel calculation as 
for nRCT, however, since the scattering is nearly adiabatic, we assume that the initial state is a 
scattering solution of the A 2 S + potential and the final state is either a bound or scattering state 
of the X 2 S + potential. Thus, the total radiative rate coefficient for an initial Ca+Yb + energy- 
normalized state | A, elm) with collision energy e, relative angular momentum quantum number I, 
and projection m is 

K elm = ^Y.\^)J(XJ\d\h,elm)\ 2 (1) 

where a is the fine- structure constant, fi is the reduced mass, and all quantities are expressed 
in atomic units. The sum is over final states /, which denote either a continuum wavefunction 
|X, e'l'm') leading to RCT or a ro-vibrational state |X, v£'m') of the X 2 E + potential leading to RA. 
The quantity uif^ is the frequency difference between the initial and final states and the operator 
d is an abbreviation for d(R)Ci g (R), where d(R) is the .R-dependent dipole moment between the 
X and A state, C\ q (R) is a spherical harmonic, R the orientation of the molecule in a space-fixed 
coordinate system, and q = m — m! gives the polarization of the emitted photon. The selection 
rules of the dipole moment ensure that \£ — £'\ < 1 and \m — m'\ < 1. 
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